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ABSTRACT 

We quantitatively examine the significance of star formation triggered in the swept- up sheU around an 
expanding H II region. If the swept-up molecular gas is sufficiently massive, new OB stars massive enough 
to repeat the triggering process will form in the shell. We determine the lower limit (Mthr) for the mass of 
the star that sweeps up the molecular gas, where at least one new star with mass Af, > Mthr forms after 
the shell fragmentation. To calculate the threshold stellar mass, Mthr, we examine how massive molecular 
shells can form around various central stars, by performing detailed numerical radiation hydrodynamics 
calculations. The mass of the photodissociated gas is generally larger than the mass of the photoionized 
gas. However, the swept-up molecular mass exceeds the photodissociated mass with a higher-mass star of 
M* > 20 Mq. The accumulated molecular mass generally increases with the stellar mass, and amounts 
to 10'*"^ M© for M* > 20 Mq with an ambient density of n ~ 10^ cm~'^. The threshold stellar mass is 
Mthr ~ 18 Mq with the star- formation efficiency of e ~ 0.1 and n ^ lif' cm^^. We examine the generality 
of this mode of run-away triggering for different sets of parameters, and found that Mthr ~ 15 — 20 Mq 
in various situations. If the ambient density is too high or the star-formation efficiency is too low, the 
triggering is not run- away, but a single event. 

Subject headings: Circumstellar matter - H II regions - ISM: molecules - STARS : formation 



1. INTRODUCTION 

In galaxies, most of the stars are born in localized ac- 
tive star-forming regions or starburst regions. Clustered 
star formation is the dominant mode in these regions (e.g., 
Evans 1999; Lada & Lada 2003), and several star clus- 
ters, including OB stars, form over a few Myrs from the 
parental giant molecular cloud (Williams & McKee 1997). 
The newly born massive stars promptly emit the UV 
{hv > 13.6 eV) and far-UV (FUV ]lleY <hv< 13.6 eV) 
radiation, and the radiative feedback by these photons oc- 
curs much earlier than the supernova explosions. The UV 
and FUV radiation from massive stars has two competing 
effects on the parental molecular cloud. One is the negative 
feedback effect on the star formation activity. The parental 
molecular cloud is ionized and heated up by the UV ra- 
diation, and the nearby star formation is quenched (e.g., 
Whitworth 1979; Franco et al. 1994; Wilhams & McKee 
1997; Matzner 2002). Diaz-Miller et al. (1998) have shown 
that the photodissociation by FUV photons is more signif- 
icant than the photoionization by UV photons. The other 
is the positive feedback effect. The next star formation is 
triggered in the compressed dense layer around the H II re- 
gion (collect and collapse scenario ; e.g., Elmegreen & Lada 
1977; Elmegreen 1989). Recent observations have provided 
dramatic snapshots of this triggering process (e.g., Dehar- 
veng et al. 2003, 2005; Zavagno et al. 2005). Keto, Ho & 
Lo (2005) have observed details of the starburst regions, 
and suggested that the propagating star formation is im- 
portant on the scale of the clouds. The character of star 
formation depends on which feedback effect dominates the 
other in the cloud. 

In our previous papers, Hosokawa & Inutsuka (2005, 
2006) (hereafter Papers I and II), we have studied the 
time evolution of the H II region, photodissociation re- 



gion (PDR), and the swept-up shell, performing detailed 
numerical calculations. We have shown that the molecu- 
lar gas is accumulated around the H II region, and that 
the gravitational fragmentation of the shell is expected 
in many cases with the homogeneous molecular ambient 
medium. This suggests that the expanding H II region 
is an efficient trigger of the star formation in the cloud. 
In this letter, we examine the net feedback effect of the 
expanding H II region and PDR around a massive star. 

2. A STARBURST MECHANISM 

2.1. Triggering Threshold Condition 

First, we focus on the positive feedback process for trig- 
gering the star formation in the shell. Let us consider 
a situation where a dense molecular shell, whose mass is 
Msh,m, forms around a massive star, and the subsequent 
star formation is triggered after fragmentation of the shell. 
The total number of newly born stars is calculated as 



M,, 



(1) 



where e is the star formation efficiency (SEE) for the 
swept-up gas and Mav is the average stellar mass, which 
is defined as, 

/■CXD I poo 

Mav = j <t>{M)M dM I j <j>{M) dM, (2) 

where (t>{M) is the initial mass function (IMF). Below, we 
use the IMF by Miller & Scalo (1979). The calculated 
value of Afav is about 0.6 Mq. The mass of the swept-up 
molecular gas, Msh,m, generally depends on the ambient 
number density and mass of the central star (see Paper 
II). If Msh.m is sufficiently large, another massive star as 
well as lower-mass stars will form in the shell. Another 
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H II region expands around this newly born massive star, 
and the triggering process wiU repeat. Furthermore, if the 
number of newly-born stars is larger than that of the pre- 
vious generation, this process causes run-away triggering 
or a burst of star formation. 

In order to examine the efficiency of this mode of trigger- 
ing, we presume that the run-away triggering begins with a 
star more massive than the threshold stellar mass, Mthr- A 
star of Mthr produces a shell that barely generates just one 
star more massive than Mthr- If the triggering cycle begins 
with a star of M* > Mthr, multiple stars of M* > Mthr 
form as "second generation" stars. The number of new- 
born stars increases as the triggering cycle advances, and 
the burst of star formation occurs. We represent the con- 
dition for the run-away triggering with simple equations. 
The number of M* > Mthr stars born in the shell is. 



-^thr = fthrN* 



(3) 



where /thr is the number ratio of M* > Mthr stars calcu- 
lated as, 



thr 



POO 



(l>{M) dM 



(j){M) dM. (4) 



We can determine Mthr from the triggering threshold con- 
dition; iVthr = 1- Combining this condition and equations 
(1) and (3), we obtain 



M 



sh,m 



M^ 



thr 



(5) 



Evaluating how Msh,m depends on the stellar mass, we can 
calculate Mthr for a given ambient number density and 
SFE by equation (5). If the molecular mass of the shell 
is smaller than Mav/e/thr, triggering is not run-away, but 
only a single event. We examine Msh,m with various stellar 
masses and ambient number densities based on numerical 
calculations. 
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2.2. Mass of the Molecular Shell v.s. Mass of the Central 

Star 

In order to evaluate the threshold stellar mass, Mthr 
for the burst of star formation, we examine how the shell 
mass, Msh,m changes with the mass of the massive star 
and the ambient molecular gas density. Even if the shell 
forms around the H II region, the accumulated gas is ex- 
posed to FUV radiation from the central star. Unless the 
FUV radiation is blocked, almost no cold molecular gas 
accumulates in the shell. Therefore, we should carefully 
estimate the mass of the cold molecular gas in the shell, 
which is available for the next star-formation episode. For 
this purpose, we have calculated the time evolution of the 
H II region, FDR, and the swept-up shell around various 
central stars in the homogeneous molecular medium (see 
Paper II for the details). In this subsection, we focus on 
the fiducial case with nn = 100 cm""^ (riHa = 50 cm~^). If 
the FUV radiation is efficiently shielded in the shell, the 
accumulated molecular mass, Msh,m increases as the H II 
region expands. However, Msh,m is actually limited by the 
following two factors: 



Fig. 1. — Top panel : Stellar mass v.s. the termination time 
of the shell in the ambient medium of nji = 100 cm~^ with the 
dissolving velocity, v^i^ = 2 km/s. The solid and broken lines rep- 
resent the dissolving time of the shell, t^i^ and the main-sequence 
lifetime of the central star, tMS respectively. The termination time 
of the shell is the smaller one between t^i-^ and tMS ■ The dotted line 
means the time when the gravitationally unstable region appears in 
the shell. The dots and open circles indicate the calculated values 
in each model. Middle panel : The mass of molecular hydrogen 
accumulated in the shell (solid line), neutral hydrogen in the PDR 
(dot-dashed line), and ionized hydrogen in the H II region (dot- 
ted line) for various central stars at the termination time. Bottom 
panel : The relative fraction of three masses presented in the mid- 
dle panel. We plot M X M^^'^ oc Mtj>{M^), where </i(M,) is the 
IMF given by Miller & Scalo (1979). In each panel, a positive feed- 
back effect is expected in the unshaded mass range, while negative 
is expected in the shaded range. 



First, the expanding shell finally dissolves due to the tur- 
bulent motion of the ambient medium. Ubiquitous su- 
personic turbulent motions in molecular clouds have been 
recognized by the observed broad emission line width. The 
corresponding velocity dispersion increases with the length 
scale as cr oc and u ~ 2 km/s at / ~ 10 pc. We de- 
fine the dissolving velocity, Udiv ~ several x km/s, and 
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the shell dissolves at the time, idiv, when the shell veloc- 
ity becomes less than Vdiv Second, OB stars evolve to 
the Wolf-Rayet phase after the main-sequence phase, and 
finally explode as supernovae. Although the Wolf-Rayet 
winds and supernova explosions must cause the dynami- 
cal expansion of bubbles, we do not consider the evolution 
after the main-sequence lifetime, tus- We evaluate the 
mass of the molecular gas in the shell at the termination 
time of the shell, min(fdiv,iMs)- The top panel of Fig- 
ure 1 presents t^w and tMS with various central stars for 
Tin = 100 cm~^ and Vdw = 2 km/s. As this panel shows, 
the dissolving time is shorter for the lower-mass star. This 
is because the dynamical time, idyn is shorter for the loss 

1 /3 

luminous central star as tdyn = -Rst/Cii c< S^y (Paper II), 
where 5uv is the UV photon number luminosity. With 
shorter tdyn, the expansion quickly decelerates and termi- 
nates earlier. The main-sequence lifetime, ^ms is shorter 
for the higher-mass star, and the expansion is limited by 
tus for > 30 M0. 

The middle panel of Figure 1 presents the mass of the 
swcpt-up molcciilar gas at the termination time. This 
panel shows that Msh,m increases with the stellar mass, 
which verifies the triggering threshold condition. The 
molecular mass amounts to Mgh 10^-5 Mq for the 
massive central stars of > 20 Mq. In these cases, the 
H II region expands to 10 — 20 pc by the termination of 
the shell. With the central star of M* < 20 Mq, the ac- 
cumulated molecular mass is significantly small. The shell 
dissolves before the accumulation of the molecular gas with 
the lower-mass star. To clarify the dominant feedback ef- 
fect, we also plot the masses of the photodissociatcd and 
photoionized gas around the star at the same time. With 
the higher-mass stars of M, > 20 Mq, the accumulated 
molecular mass is much larger than both the ionized and 
photodissociatcd masses. Therefore, such massive stars 
give have a positive feedback effect on the clouds. With 
the lower-mass stars of M* < 20 Mq, on the other hand, 
the photodissociatcd mass is much larger than the ionized 
and accumulated molecular masses. In these cases, the net 
feedback effect is negative, and the photodissociation by 
FUV photons is the dominant negative process, as argued 
by Diaz-Miller et al. (1998). The bottom panel of Figure 1 
presents the masses multiplied by the population of excit- 
ing stars. When the higher-mass stars of Af, > 20 Mq are 
formed in the cloud, the positive feedback effect can domi- 
nate, despite a large number of lower-mass stars. Whether 
the net feedback effect is negative or positive depends on 
the mass of the most massive star in the cloud; positive for 
M, > 20 Mq, and vice versa in the current fiducial case. 



2.3. Threshold Stellar Mass for Starburst 

Using the calculated Msh,m for various central stars, we 
determine the threshold stellar mass, Afthr for the burst 
of star formation. For example, in the fiducial case with 
nil = 100 cm~^ and Vdw = 2 km/s, the threshold mass 
calculated by equation (4) is Mthr ^ 18 Mq with e = 0.1. 
This means that only a star more massive than 18 Mq can 
generate a shell massive enough so as to produce at least 
one star as massive as the original star. 
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Fig. 2. — Threshold stellar masses, M^^^ for the run-away trig- 
gering with different parameters. The horizontal axis is the ambi- 
ent number density, and the different symbols denote different sets 
of dissolving velocity, t'^iv and star formation efficiency, c. Some 
symbols disappear at higher ambient densities, which means that 
the triggering threshold condition (eq.(5)) is not satisfied with the 
corresponding set of pcirameters. In these cases, triggering is not 
run-away, but a single event. 

Note that the threshold mass is comparable to the stellar 
mass over which the character of feedback changes from 
negative to positive. In this case with M* > 20 Mq, there- 
fore, the net feedback effect is positive, and strong enough 
to cause the run-away triggering. 

We have examined the possibility of run-away trigger- 
ing by calculating Mthr with different parameters (Fig. 2). 
For example, Mthr increases with decreasing SFE. This 
is because the more massive molecular shell is needed to 
compensate for the lower SFE. With e = 0.03 in the fidu- 
cial case, we obtain Mthr ~ 21 Mq. If the SFE is as low 
as e < 0.01, however, equation (4) is not satisfied with any 
Mthr- This is because the increase in Msh,m grows satu- 
rated at about 10^ Mq as the mass of the central star 
increases (sec Fig.l). Although the net feedback effect is 
positive with M* > 20 Mq, triggering is not run-away, but 
only a single event. 

The threshold stellar mass also increases with the higher 
ambient density. In the denser ambient medimn, the size 

of the H II region shrinks as Rg^ ocn^ ' . The total swept- 
up mass at the termination time also becomes smaller fol- 
lowing the scaling relation, Mgh oc n^^, with the same 
central star. With the higher ambient density, therefore, 
Mthr increases to produce the massive shell which enables 
the run-away triggering. Since the increase in Msh,m be- 
comes saturated for the higher-mass stars, the run-away 
triggering becomes impossible unless the SFE is very high 
in the denser ambient medium. With uh = 10^ cm~'^ and 
Wdiv = 2 km/s, for example, an SFE as high as e ~ 0.3 is 
needed for the run-away triggering. 

Finally, the higher dissolving velocity reduces the signif- 
icance of triggering. The expanding shell dissolves earlier 
with higher Vd\v, and the H II region sweeps up only the 
smaller region. Furthermore, the molecular fraction in the 
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shell is preferentially reduced with the higher Wdiv During 
the dynamical expansion of the H II region, the column 
density of the shell, Nsh increases as Ngh ~ nu RiF{t)/5. 
Therefore, the smaller H II region has a shell with lower 
column density, though the incident FUV radiation is 
strong owing to the poor geometrical dilution. Since 
the FUV radiation shielding is inefficient, it is hard for 
molecules to accumulate in the shell. With Wdiv = 4 km/s 
and e = 0.3, for example, the run-away triggering is possi- 
ble only with the low ambient density of nu ^ 10^'^ cm~^. 

3. DISCUSSIONS AND CONCLUSIONS 

In this letter, we suggest a possible mode of self- 
propagating star formation in molecular clouds, which 
should be examined in further studies. Although we have 
adopted the homogeneous ambient density, for example, 
real clouds show clumpy structure with a turbulent ve- 
locity field. Recently, Mellcma ct al. (2006) have calcu- 
lated the dynamical expansion of the H II region in turbu- 
lent molecular clouds. They have shown that the clumpy 
medium leads to an irregular ionization front. The ra- 
dial position of the ionization front approximately agrees 
with that in the homogeneous medium on average with the 
mean density, but significantly depends on the radial direc- 
tion. More theoretical studies and multi-dimensional nu- 
merical simulations of the deformation and fragmentation 
processes of the swept-up shell are needed. Observational 
estimation of the SFE in the shell might be interesting and 
useful for the present model. 

The stellar wind from the massive star is another omis- 
sion in our evaluation. The expanding wind-driven bubble 
can modify the dynamics of the ISM around massive stars. 
When the bubble pressure is much higher than the ambi- 
ent pressure, the time evolution of the bubble size and 
pressure is given by, 

''b(«) = (j^)"'"iyv"v'». (6) 

where Lw is the wind mechanical luminosity (Weaver et 
al. 1977). Here, we evaluate the significance of the wind- 
driven bubble by comparing the ratio between the bubble 
pressure and H II pressure (Pn = PoC'n) at i?b = -Rst- 
With equations (6) and (7), 
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where we define /l = L^/L^t and Lst = 47ri?g^poC'n (Mc- 
Kee et al. 1984). If /p < 1, the bubble is confined by the 
ambient H II pressure before reaching the initial Stromgren 
radius. The numerical value of /l can be estimated by 
/l = 0.8{Suv/10^^s-^f/'^nll^ (Abbott 1982). With the 
number density of nn = 100 cm~^, /p ~ 0.4 for a 20 Mq 
star and /p ~ 1.5 for a 50 M© star. Therefore, the stellar 
wind only slightly increases the initial pressure of the H II 
region with the star of M* < 50 Mq. The effect of the 
wind-driven bubble becomes significant with a very mas- 
sive star (~ 100 Mq), highly clustered massive stars, or in 
higher density media. When a wind-driven bubble domi- 
nates the dynamics, the bubble sweeps up the larger region 
at higher velocity than the expansion only due to the H II 
region overpressure. This situation may even promote the 
positive feedback effect, and should be separately studied 
in detail. 

Finally, we summarize the results. First, we have formu- 
lated the conditions for the positive feedback process by 
which the triggered star formation continues in the swept- 
up shell. The threshold stellar mass, Afthr, is defined as 
the mass of the central star that sweeps up the molecular 
gas in the shell, where at least one new star with mass 
M* = Mthr forms after the shell fragmentation. In order 
to evaluate Mthr, we have calculated how massive molec- 
ular shells can form around various central stars. In the 
fiducial case with nn = 100 cm~'^ and v^w = 2 km/s, the 
swept-up molecular mass increases with the stellar mass, 
and amounts to Mgh^m ^ 10''"^ Mq with the massive star 
of M„ > 20 Mq. We have also calculated the photoionized 
and photodissociated masses to clarify the net feedback ef- 
fect. The negative effect of photodissociation is important 
with the lower-mass star of M* < 20 Mq, but is dominated 
by the positive effect with the higher-mass star. We have 
calculated the threshold mass as Mthr ~ 18 Mq in the 
fiducial case with e = 0.1. We have also calculated Mthr 
with different parameters, and examined the generality of 
run-away triggering. We have found Mthr ~ 15 — 20 M© in 
various situations. The triggering process is only a single 
event with the higher ambient density, higher dissolving 
velocity, or lower SFE. 
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